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ABSTRACT 

We present time- resolved optical spectroscopy of the dwarf nova GW Librae during 
its rare April 2007 super-outburst and compare these with quiescent epochs. The data 
provide the first opportunity to track the evolution of the principal spectral features. 
In the early stages of the outburst, the optically thick disc dominates the optical and 
the line components show clear orbital radial velocity excursions. In the course of 
several weeks, optically thin regions become more prominent as strong emission lines 
replace the broad disc absorption. 

Post-outburst spectroscopy covering the I-band illustrates the advantages of Ca ii 
relative to the commonly used Balmer lines when attempting to constrain binary 
parameters. Due to the lower ionisation energy combined with smaller thermal and 
shear broadening of these lines, a sharp emission component is seen to be moving 
in between the accretion disc peaks in the Ca ii line. No such component is visible 
in the Balmer lines. We interpret this as an emission component originating on the 
hitherto unseen mass donor star. This emission component has a mean velocity of 
^ —15 ±5 km s~^ which is associated with the systemic velocity 7, and a velocity semi- 
amplitude of Kem = 82.2 ± 4.9 km s~^. Doppler tomography reveals an asymmetric 
accretion disc, with the S-wave mapping to a sharp spot in the tomogram with a 
velocity consistent to what is obtained with line profile fitting. A centre of symmetry 
analysis of the disc component suggests a very small value for the WD orbital velocity 
Ki as is also inferred from double Gaussian fits to the spectral lines. 

While our conservative dynamical limits place a hard upper limit on the binary 
mass ratio of q < 0.23, we favour a significantly lower value near q ^ 0.06. Pulsation 
modeling suggests a WD mass ~ IMq. This, paired with a low mass donor, near 
the empirical sequence of an evolved CV close to the period bounce, appears to be 
consistent with all the observational constraints to date. 

Key words: binaries: spectroscopic - novae, cataclysmic variables - stars: individual: 
GW Lib 



1 INTRODUCTION 

Cataclysmic variables (CVs) are semi-detached interacting 
binary systems containing a white dwarf (WD) primary and 
a late-type main-sequence secondary (MS) which transfers 
mass on to the primary through Roche Lobe overflow, see 
I Warned (1995) for a review. Strong and broad emission lines 
are a common observational feature in the majority of CVs, 
probing the dissipative processes within the accretion flow. 
The hydrogen Balmer series usually dominate the spectrum, 
and in the vast majority of spectroscopic studies of CVs 
serve as the main proxy for determining the system param- 
eters. However, relatively unexplored spectral features re- 
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main, such as the Ca 11 triplet lines in the I-band. Common 
emission line sources include the accretion disc, the interac- 
tion point of the gas stream with the disc (hot spot) and the 
irradiated face of the secondary star. 

GW Lib is a dwarf n ova (DN) which was discovered by 
IConzalez fc Maze! l|l983h as it went onto its first recorded 
outburst and rapidly increased in brightness by 8 magni- 
tudes. It remained in it s quiescent sta t e afte r that period 
of activity until 2007. IVan Zvl et ai] |20oO) found three 
periodicities in the quiescent photometric light curve at 
236, 376 and 650 s, making the primary in GW Lib the 
first multi- period non-radial pulsati ng WD in an accret- 
ing binary. [Thorstensen et alT l|2002l ) measured the radial 
velocity curve using Ha and determined a rather short 
orbital period of 76.79 ± 0.02 min. In the latest Ritter 
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& Kolb catalogue, only 8 systems with a smaller period 
are known when exclud ins the AM CVn s ystems which 
have H-deficient donors (jRitter fc Kolbir2003h . implying an 
evolved binary close to the period minimum. From model 
fits to the WD absorption troughs and mean spectra, var- 
ious temperatures have been found. UV sp ectra suggest 
Two = 14, 700 K with a WD mass of O.6M0 (|Szkodv et all 
bOOd ) or a dual temperature WD with Tiow = 13, 300 K and 
Thigh = 17, 000 K. The latter model can only explain the 
observed UV/optical pulse am plitude ratio in GW Lib (see 
table 2 in lSzkodv et~aLll2002l ) as compared with single pul- 
sating ZZ-C eti stars, if the WD mas s of GW Lib is larger 
thanO.SMr^. lThorstensen et all l|2002l ) fit the optical Balmer 
line profiles to find T = 13, 220 K and log g — 7.4. However, 
all temperature estimates place the WD outside the ZZ Ceti 
instability strip for single WD pulsators. A mass determina- 
tion based on the detected pulsation periods combined with 
distance and UV flux constra ints suggests a WD mass of 
1.02Mq l|Townslev et al.ll2004h . 

On 2007 April 12, amateur astrono mers reported the 
sudden and rapid brightening of GW Lib (|Templetonll2007l ') 
indicating renewed outburst activity, 24 years after the dis- 
covery outburst. The outburst had an amplitude of ~ 9 mag- 
nitudes and lasted for several weeks, see Figure [T] 

We present the first time-resolved optical spectroscopy 
of GW Lib in (super)outburst (Section [Sj alongside time- 
resolved optical spectroscopy during quiescence before and 
after the outburst (Section |4]). We present the possible sys- 
tem parameters in Section [5] and the discussion in Section 

El 



2 OBSERVATIONS AND REDUCTION 

We briefiy describe the facilities used below. Full details of 
the observations can be found in Table [T] 

2.1 Telescopes and Instruments 

The telescopes used were the Magellan Telescopes oper- 
ated by the Carnegie Institution of Washington at the Las 
Campanas Observatory in Chile, the Smithsonian Astro- 
physical Observatory's Tillinghast Telescope at the Fred 
Lawrence Whipple Observatory, located on Mount Hopkins 
near Amado, Arizona and the William Herschel Telescope 
operated by the Isaac Newton Group on the island of La 
Palma. 

2.1.1 Magellan Telescope 

Time-resolved optical spectroscopy was acquired during qui- 
escence with the Baade 6.5-m telescope equipped with the 
Boiler and Chivens Spectrograph (b&c) on 2004 June 6 and 
7. The B&C was used with a Marconi 2048 x 515 CCD with 
a 13.5 microns pixel size and a 1200 lines mm"'^ grating 
covering the spectral interval 3500 — 5100A. A slit width of 
0.8 arcsec gave a dispersion of 0.79A pixeP^ and a resolution 
of 2.0A. 

On 2007 June 19 the telescope was equipped with the 
Inamori-Magellan Areal C amera and Spectrograph (imacs: 
iBigelow fc Dressier! |2003| ) to acquire time-resolved spec- 
troscopy of GW Lib. Imacs was used with the f/4 camera 



and the long-slit-mask with a 0.7 arcsec slit width and a 600 
lines mm~^ grating centred at 5550A. The spectra were dis- 
persed along the short-axis of four SITe CCDs in the IMACS 
mosaic detector. The CCD detectors were binned two by two 
during the observations. This instrumental setup provided 
a dispersion of 0.76A pixel"^ and a spectral resolution of 
1.29Afwhm in the spectral interval 3950 - 7100A. 

2.1.2 Tillinghast Telescope 

The 1.5-m refiector Tilli nghast Telescope (tt ) equipped 
with the FAST instrument l|Fabricant et al.ll 19981 ) acquired 7 
epochs of time-resolved spectroscopy and 16 epochs of single 
spectroscopic frames in service mode. For the time-resolved 
observations, the FAST instrument was equipped with a 2688 
X 512 UA STA520A CCD chip with a 15 microns pixel size. 
A grating of 1200 lines mm~^ covering 4150 — 5200A and 
a 1.5 arcsec slit width was used to deliver a dispersion of 
0.38A pixel"'^ and spectral resolution of 0.86A fwhm. 

For the single frames the telescope was used in service 
mode with a grating of 300 lines mm^^ and an aperture 
of 3.0 arcsec to cover the spectral interval 3300 — 7600A. 
The pixels were binned two by two to give a dispersion of 
I.47A pixel"'^ and a resolution of 5.9A. 

2.1.3 William Herschel Telescope 

On two successive nights, three months after the 2007 out- 
burst, we acquired time-resolved spectroscopic data using 
the 4.2-m William Herschel Telescope (wht) in combination 
with the two-armed Intermediate dispersion Spectrograph 
and Imaging System (isis). The blue arm of the spectro- 
graph was equipped with a 4096 x 2048 EEV12 CCD. The slit 
width used was 1.0 arcsec and used with the R1200B grat- 
ing with 1200 lines mm~^ providing a wavelength coverage 
of 4200 - 5OOOA. This setup delivers a dispersion of 0.224A 
pixel" ^. The red arm was equipped with a 4096 x 2048 
Red-|- CCD with 15 microns pixels. Using the R1200R grat- 
ing with 1200 lines mm~^ the unvingetted wavelength cov- 
erage was 8050 — 88OOA with a dispersion of 0.243A pixeP^. 
The spectral resolution for both arms was O.62A. 

2.2 Reduction 

The different sets of spectroscopic data were reduced using 
several reduction packages. 

2.2.1 FAST service mode 

The single frame spectra obtained with the FAST spectro- 
graph in service mode were extracted using the spectral 
extraction pipeline provided by the Telescope Data Center 
at the Harvard-Smithsonian Center for Astrophysics. This 
tailored pipeline is b ased on the iraf procd package, see 
iTokarz fc Rolil l| 19971 ). 

2.2.2 Magellan IMACS 

The IMACS data were bias and fiat-field corrected in the stan- 
dard way using iraf. The spectra were extracted from each 
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Table 1. Observations of GW Lib 



Telescope 


Instrument 


Date 


HJD 

(start) 


Exp time 

(s) 


n spectra 


Urbital 
coverage 

l-'obsZ-'orbJ 


A - range 
(A) 


Seeing 
( ) 


Magellan 


B&C 


06/07/2004 


2453193.491 


60 


62 


1.00 


3500 - 5100 


0.5-1.0^ 






07/07/2004 


2453194.471 


60 


66 


1.02 




0.7-1.2" 


TT 


FAST 


14/04/2007 


2454204.859 


60 


1 




3300 - 7600 


1-2 






14/04/2007 


2454204.861 


30 


2 






1-2 






20/04/2007 


2454210.859 


15 


2 












09/05/2007 


2454229.833 


120 


1 






2 






09/05/2007 


2454229.836 


300 


1 






2 






11/05/2007 


2454231.817 


120 


2 






2 






15/05/2007 


2454235.813 


480 


1 












18/05/2007 


2454238.894 


600 


1 


- 




- 






20/05/2007 


2454240.811 


480 


1 












23/05/2007 


2454243.787 


720 


1 


- 




- 






09/06/2007 


2454260.739 


720 


1 












16/06/2007 


2454267.774 


1200 


1 












24/06/2007 


2454275.691 


1200 


1 


- 




- 


TT 


FAST 


14/04/2007 


2454204.886 


200 


34 


1.86 


4150 - 5200 


1-2 






15/04/2007 


2454205.861 


60 


110 


1.75 




1-2 






20/04/2007 


2454210.882 


60 


70 


1.46 




1-2 






23/04/2007 


2454213.858 


120 


37 


1.68 




1-2 






10/05/2007 


2454230.833 


300 


18 


1.75 




2 






12/05/2007 


2454232.833 


600 


2 


0.16 




2 






17/05/2007 


2454237.784 


300 


22 


1.94 




1-2 


Magellan 


IMACS 


18/06/2007 


2454270.531 


30 


73 


1.09 


3950 - 7100 


0.6 


WHT 


ISIS 


24/07/2007 


2454306.365 


240 


19 


1.08 


4200 - 5000/8050 - 8800 


1 






25/07/2007 


2454307.350 


240 


25 


1.35 




1-2 
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CCD frame with the iraf kpnoslit package. The pixel-to- 
wavelength calibration was derived from cubic spline fits to 
HeNeAr arcs acquired during the observations. Fitting 10 - 
20 lines per frame gave a RMS < O.OI7A per chip. 

2.2.3 All others 

All other data sets were reduced in the following way. The 
average bias and flat-field correction was carried out using 
the FIGARO package from STARLINK and nightly average bias 
and tungsten frames. Pamela was used for the optimal ex- 
traction of the spectra (|Marshlll989l ). Regular CuAr arc lamp 
exposures allowed us to establish an accurate wavelength 
scale for each spectrum through interpolation between the 
nearest arcs in time. Each arc frame has been fitted with a 
5th order polynomial to 10-40 lines to give a typical RMS 
of 0.1 pixel. The individual spectra were normalised to the 
continuum level using a spline fit to selected continuum re- 
gions. 

3 GW LIB IN OUTBURST 

The April 2007 outburst of GW Lib triggered our target of 
opportunity program at the TT. Both time-resolved spec- 



troscopy and single frame observations of GW Lib were ob- 
tained during this period of activity. See Table [T] for speciflc 
details. 



3.1 Spectral evolution 

In Figure [5] average representative normalised spectra dur- 
ing different stages of the outburst of GW Lib are presented. 
The bottom spectrum is during quiescence 3 years before the 
2007 outburst, the second spectrum is a day after outburst 
peak, the third during slow decline, the fourth during rapid 
decline and the fifth during the final decline towards quies- 
cence. For reference, the top frame in Figure [T] shows the 
outburst light curve in the V-band, as provided by AAVSO 
and its network of observers. The dates corresponding to the 
spectra plotted in Figure [5] are marked. The pre-outburst 
quiescent spectrum is given an arbitrary negative HJD for 
plotting purposes. 

The single frame spectra at different nights through- 
out the outburst make it possible to trace the evolution of 
various spectral features. We have measured the equivalent 
widths (ew) of the lines as a function of the outburst. For the 
Balmer lines, the overall continuum is normalised but any 
underlying absorption within the line was not removed. This 



4 L. van Spaandonk et al. 



i- ' ' ' 

1 

fl 

II II 




^ : 

; 










*^ 24 07 2007 -5 

^f^l>)lf■l^^^<^f^,l^4.l^l^|^)trtl^^l^l*»%i> 






12 05 2007 -4 






23 04 2007 -3 




A ^ 


15 04 2007 - 2 




——J ■—' • ' : 07 07 2004 - 1 



Wavelength 

Figure 2. Comparison of average representative normalised spectra of GW Lib sampling the spectral evolution before and during 
the 2007 super-outburst. The April 2007 spectra have been multiplied by a factor of 4 to magnify the spectral features and allow 
for a better comparison. The vertical dotted lines annotate the ±2500 km s~^ window for the Balmer lines and the ±1500 km s~^ 
band for helium lines while the horizontal dotted lines the normalised continuum flux. 




Figure 1. Time evolution of the equivalent widths of the prin- 
cipal spectral lines. The July 2004 quiescent epoch is plotted at 
an arbitrary negative HJD for comparison. In the bottom frame: 
solid line for H/3 and dashed line for H7. In the middle panel: 
solid line for He I 492lA, represents the general evolution of the 
He I lines and dashed line for He II 4685A. The top panel shows 
the V magnitude of the outburst, as provided by the AASVO where 
the numbers indicate position of the spectra from Figure [2] which 
are discussed epochs in Section [XT] 



results in positive (ew) values during absorption-dominated 
epochs and negative values otherwise. For consistency, the 
window for these lines is ±2500 km s~^ for all epochs, 
which is limited by the proximity of the helium lines. For 
the He lines, a slightly narrower window of ±1500 km s~^ 
gives a consistent and comparable EW measurement, avoid- 
ing neighbouring lines. The spectra are locally normalised 
to avoid any contamination from underlying absorption due 
to nearby Balmer lines. Both regions are marked in Figure 
[2] together with the continuum level. The bottom panel in 



Figure [T] shows the EW evolution for the Balmer lines and 
the middle panel the same for the helium lines. The actual 
values for all lines and at all observed epochs are given in Ta- 
ble [2] We note that the spectral lines display rather complex 
changes as different components come and go. A single (ew) 
value clearly does not fully capture these details. Therefore, 
to complement the measured (ew), a more detailed anal- 
ysis of the separate components within the line profiles is 
provided below. 

2004 July 7-1 During deep quiescence, 3 years before 
outburst, the average spectrum of GW Lib was typical for 
a CV spectrum with a low accretion rate. The accretion 
disc has very narrow Balmer emission lines which show a 
dip in their centres refiecting an underlying double-peaked 
profile (not visible on the scale of Figure [2]). This profile 
shape is expected for optically thin accretion discs viewed 
at a low inclination angle i. The H/3 profile is fitted with two 
Gaussians simultaneously to the absorption trough and the 
narrow emission line component. The latter gives a fwhm 
of 7.82 ± 0.56A corresponding to ~ 482 km s^^ (all er- 
rors quoted in this paper are 3cr errors) . Absorption troughs 
flank the emission features with a FWZI of ^ 20, 000 km s^^, 
and can only originate from the WD as these velocities are 
too high for Kepler velocities in the accretion disc. This in- 
dicates a small accretion luminosity and a very low mass 
transfer rate such that the primary WD is visible. The mea- 
sured EW of the Balmer lines are dominated by the WD 
absorption. The He I lines are all in emission. A small, but 
significant, amount of He 11 is present in emission with an 
EWof-0.34±0.18A in the average spectrum. Since the disc 
and companion star are expected to be too cool during qui- 
escence to produce much He 11 emission, this line emission 
may be due to the relatively hot WD ionising gas near it. 
The feature is unfortunately too weak to permit any time- 
resolved analysis. We also note that an absorption line at 
Mg 11 4481. 15A is present with an EW of 0.24±0.09A mea- 
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sured in a ±300 km window and a fwhm of 3.18±0.35A. 
This is likely formed in the photosphere of the WD and 
such metal lines may be expected in accreting systems were 
metal-rich gas is deposited onto the WD surface. The deter- 
mination of the gravitational redshift of such lines could in 
principle provide a direct measurement of the mass of the 
WD. However, our resolution prohibits a reliable measure- 
ment of the velocity of this weak feature as it is blended 
with the nearby He l line. 

2007 April 15-2 Two days after the rise into outburst 
started, the average spectrum of GW Lib shows broad ab- 
sorption troughs in the Balmer lines. Given the 9 magnitude 
increase due to accretion these absorption features can no 
longer be associated with the WD but suggest an origin in 
the now optically thick accretion disc, supported by a fwhm 
of 23.05 ±0.36A. The lines also show strong emission peaks, 
see spectrum 2 in Figure [2 This is al s o see n in the Ha 



emission profile obtained bv lHiroi et al.l (|2009l ) around this 
time. With a fwhm of only 3.14 ± 0.66A, corresponding to 
~ 194 km s~^, the peak in H/3 has less than half the width 
of the emission during quiescence. The He I lines are all in 
absorption, with a hint of a very narrow emission peak at 
the centre, while the He ll line at 4685A is in emission at 
EW = -0.25 ± 0.19A. Finally Fe is marginally detected in 
absorption. We will revisit the discussion on the origin of 
these lines in section [321 



2007 April 23-3 Ten days after the beginning of the 
outburst, GW Lib slowly started to fade and the spectrum 
evolved. The strong and narrow emission in the H/3 and 
H7 lines decreased in strength, increasing the overall EW 
slightly. The EW of the He I lines increases similarly suggest- 
ing the same narrow emission was indeed present. The accre- 
tion disc is still visible as an optically thick disc in absorption 
in the Balmer lines with a fwhm of 30.46 ±0.26A. The He n 
emission decreases slightly to an EW of —0.32 ± 0.06A. 

2007 May 12 - 4 A month after outburst, as the system 
dropped from the outburst plateau (Figure [T]) and the lu- 
minosity decreased rapidly, the spectrum had changed dras- 
tically, see spectrum 4 in Figure [S] The accretion disc con- 
tribution to the lines changed from an optically thick flow 
back into a chiefiy optically thin disc, producing shallow 
double-peaked emission profiles in the Balmer line without 
the high velocity absorption, leading to a negative EW. The 
emission component in H/3 is visible again and now has a 
FWHM of 8.59±0.03A, close to the value of the pre-outburst 
emission line. As can be seen in Figure [T] the EWs of all 
lines make a small jump. WZ Sge type DNe often show brief 
brightenings short ly after their main o utburst, referred to 
as echo-outbursts (|Patterson et al.ll2002l ') . It remains unclear 
due to the sparse light curve sampling whether the observed 
behaviour in the lines is related to similar events. 

2007 June 19 The Magellan spectra from June 2007 re- 
semble the WHT July spectra in shape, see top spectrum in 
Figure[2]for the average spectrum. In addition, these spectra 
also cover the Ha and the He i 5876. 62A line. An interesting 
feature next to the strong He I line is the Sodium doublet 
(Na D 5889.95A-5895.92A) seen in emission, see Figure El 




6540 6550 6560 6570 6580 5860 5870 5880 5890 5900 5910 
Wavelength ' Wavelength 



Figure 3. Normalised emission line profiles of Ha (left) and the 
He I - Sodium doublet (right). The solid line represents in both 
cases the average profile over 1 orbital period. The dashed lines 
are the various Gaussian profiles fits and the dotted line is the 
sum of all single Gaussians. As is clear from these fits, the Sodium 
doublet has, most likely, a very similar origin as the Ha and the 
He I lines. 



Due to the close proximity of GW Lib, interstellar absorp- 
tion is expected to be weak, and indeed no evidence for dif- 
fuse interstellar band (dib) absorption features can be found 
in the spectrum. Instead, the doublet is dominated by broad 
emission components. To strengthen this case, the Sodium 
doublet is also present in emission at 6154.23A-6160.76A. In 
the orbital average spectrum, the feature can be fitted with 
two double Gaussians (plotted in Figure [3} with the peak 
separation set to 250 km s~^. The Ha line has the same 
separation where for the nearby He I the line separation is 
350 km s~^. The offset of the centre of the double-peaked 
profiles from the rest wavelength for all lines is, within er- 
rors, — 10 ± 5 km s~^. These similar fitting values suggest a 
shared origin in the accretion fiow for all lines. Time-resolved 
analyses of the He l and the Balmer lines give similar values 
compared to those in the July 2007 epoch (see section 14.21 
and Table ID). 

2007 July 24-5 Three months after outburst, during 
the slow decline towards quiescence, the system was still 1- 
2 mags brighter than the pre-outburst magnitude, and the 
double-peaked profile from the disc proved stronger against 
the dropping continuum compared to the previous phase 
(compare spectrum 4 to spectrum 5 in Figure^). The system 
has cooled sufficiently such that the WD absorption again 
flanks the emission lines in the Balmer lines, decreasing the 
EW of H7. Since both the disc and the WD are expected 
to be slightly hotter and brighter compared to the pre- 
outburst quiescence state it does not yet resemble spectrum 
1 in Figure O For comparison, the WD in WZ Sge was still 
cooling after 17 m onths of post-outburst observations, see 
iLong et al.l (|2004l ). The He I lines are all in emission and the 
He II emission line has an EW of —0.52 ± 0.24A. The fwhm 
of the emission component in H/3 is 7.65 ± 0.03A and has 
thus reached the value observed in the pre-outburst epoch 
even though the emission line flux itself is still much larger. 



3.2 Time-resolved spectra 

Several epochs (see Table [T| of time-resolved spectroscopy 
allows us to characterise the line behaviour in more de- 
tail. During outburst, these time-resolved spectra showed 
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Table 2. Equivalent widths in A ±3(T error, of various emission lines as function of the outburst. The Balmer lines are measured 
in a ±2500 km s^^ window, where the helium lines are measured in a ±1500 km s~^ window 



Epoch 


HJD 


H/3 


H7 


Hel 4387 


Hel 4471 


Hell 4685 


Hel 4921 




(start) 




















1 


240oly4.4 ( i 


y.bi ± U.lo 


1 HQ _L n 1 


n on _L n 1 


n '7/1 _L n 1 
-U. r4 ± u.lz 


-0.34 


± 


0, 


,18 


-U. (O ± U.z4 


2 


2454204. 8o9 


3.32 ± 0.16 


3.55 ± 0.17 


0.04 ± 0.18 


0.47 ± 0.18 




-L 

m 


n 

u. 


1 Q 


0.10 ± 0.19 




2454210.859 


6.79 ± 0.37 


6.75 ± 0.37 


0.04 ± 0.30 


1.19 ± 0.30 


-0.35 


± 


0, 


.30 


0.29 ± 0.31 


3 


2454213.858 


7.78 ± 0.07 


7.32 ± 0.09 


0.50 ± 0.06 


1.03 ± 0.07 


-0.32 


± 


0, 


.06 


0.28 ± 0.03 




2454229.833 


-16.61 ± 1.19 


-13.37 ± 1.26 


0.29 ± 0.91 


-2.84 ± 0.97 


-0.59 


± 


0, 


,87 


-1.23 ± 0.89 




2454231.817 


-16.80 ± 1.23 


-12.85 ± 1.33 


-0.00 ± 0.98 


-1.30 ± 0.98 


-1.21 


± 


0, 


,92 


-2.15 ± 0.94 


4 


2454232.833 


-17.37 ± 0.90 


-9.77 ± 0.95 


-0.96 ± 0.62 


-0.01 ± 0.62 


-1.11 


± 


0, 


,65 


-0.98 ± 0.56 




2454235.813 


-15.07 ± 0.74 


-9.04 ± 0.78 


-0.36 ± 0.58 


-0.84 ± 0.57 


-1.33 


± 


0, 


,55 


-1.58 ± 0.56 




2454238.894 


-16.30 ± 1.01 


-9.06 ± 1.48 


-0.83 ± 1.07 


-1.06 ± 0.95 


-0.68 


± 


0, 


,80 


-2.17 ± 0.74 




2454240.811 


-19.30 ± 0.84 


-10.46 ± 0.94 


-1.15 ± 0.71 


-1.43 ± 0.66 


-0.62 


± 


0, 


,61 


-2.33 ± 0.62 




2454243.787 


-20.15 ± 0.82 


-9.43 ± 0.82 


-1.51 ± 0.64 


-1.32 ± 0.62 


-0.60 


± 


0, 


,59 


-2.96 ± 0.62 




2454260.739 


-17.25 ± 0.95 


-4.94 ± 0.93 


-1.09 ± 0.74 


-0.29 ± 0.68 


-0.82 


± 


0, 


,68 


-2.56 ± 0.73 




2454267.743 


-13.28 ± 2.05 


-2.88 ± 2.15 


-0.57 ± 1.67 


-0.74 ± 1.54 


-1.84 


± 


1, 


,53 


-3.39 ± 1.64 




2454270.531 


-15.65 ± 0.41 


-4.14 ± 0.57 


-0.65 ± 0.46 


-2.04 ± 0.42 


-1.40 


± 


0, 


,30 


-3.52 ± 0.27 




2454275.691 


-16.01 ± 1.58 


-5.62 ± 1.54 


-1.32 ± 1.28 


-3.62 ± 1.26 


-0.97 


± 


1, 


,16 


-4.03 ± 1.31 
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2454306.355 


-17.16 ± 0.35 


-6.76 ± 0.27 


-0.75 ± 0.24 


-1.47 ± 0.21 


-0.52 


± 


0, 


,24 


-2.56 ± 0.27 



considerable orbital dependence in the Balmer, He I and 
He 11 line profiles. We binned the spectra in bins of l/20th 
of the orbital period us ing P = 0.05332 ± 0.00002 days 
iThorstensen et al]l2002l '). The Balmer profiles were fitted 
with a double Gaussian with the peak, fwhm and the offset 
from the rest wavelength as free parameters. The He i and 
He 11 profiles are fitted with a single Gaussian. All fits show 
that the peak and fwhm of the line components are constant 
within 10 per cent of their average values and as these par- 
ticular parameters show no sign of significant orbital depen- 
dence, they were thus fixed to their mean values for further 
analysis, keeping their velocity offsets as free parameters. In 
addition to our profile fits, all results were checked against 
a tra ditional double Gaussian analysis l|Schneider fc Youngj 
in combination with a diagnostic diagram. 

The radial velocity curves of our line components were 
fitted with the function V{(j>) ~ ^ — K sm{27v<f) — 27r(j!)o) and 
individual velocity data points were weighted according to 
their errors as derived from the profile fits. These errors were 
scaled such that the radial velocity curve fit has a goodness 
of fit close to 1. Here 7 is the systemic velocity, K the radial 
velocity semi-amplitude of the absorption and/or emission 
line and (j>o the phase offset relative to the ephemeris. The 
ephemeris zero point will be derived later in this paper (sec- 
tion |4]) and is used throughout this paper when calculating 
orbital phases. Due to the lack of precision in the binary 
orbital period, the ephemeris cannot be extrapolated very 
far in time and thus the binary phase is arbitrary except for 
the July 2007 epoch where the zero-point is measured. 

We discuss the results of the fits for the H/3 and He 11 
profiles using the 2007 April 15 data as this set has the high- 
est signal to noise, the best time resolution and the most fea- 
tures present for analysis. The behaviour of the lines during 
the other nearby epochs were the same, albeit determined 
with lower precision. The results can be found in Figure |4] 
and Table |3] Both H/3 and H7 have similar amplitude and 
phase fits. The deep and broad absorption troughs can be 
identified with the accretion disc. However, they do not trace 
the WD movement as the velocity semi- amplitude is too 



large at 51.3±6.5 km s~ . The radial velocity of the primary 
WD is likely much smaller. Interestingly, this amplitude is 
consistent with that of the double-peaked emission profile 
during quiescence. On 2004 July 6, the velocity curve fit to 
H/3 emission gives a semi-amplitude of 45.1 ± 7.3 km s~^. 
This similarity suggests that the distortions in the disc lead- 
ing to the observed radial velocity curve of both the disc 
emission during quiescence and the broad disc absorption 
during outburst are rather similar. Unfortunately, due to 
the uncertainty in the binary ephemeris, we cannot compare 
their relative phasing. Similar amplitudes, systemic veloci- 
ties and phase zeros were found for the He I lines, suggesting 
that their profiles trace the same components. In contrast, 
the He 11 line moves in anti-phase and has half the ampli- 
tude. This highlights the difficulty with associating either 
of these radial velocity amplitudes with that of the primary 
WD. 

In both the H/3 and H7 profile, a second Gaussian is 
fitted to the narrow emission components in their cores. This 
component shows a remarkably small velocity amplitude of 
3.0 ± 1.6 km s~^ (see Table At this epoch, the arc fit is 
good to O.OO3A at the position of H/3 which is smaller than 
1/lOOth of the 0.38A/pixel scale, la errors on individual 
velocity measurements are much larger (~0.085 pixels) than 
this and are thus dominated by the statistical error. These 
errors are then propagated to provide the formal errors on 
the fit parameters as listed thus suggesting that this small 
amplitude is measured to a significant degree. 

This strong central emission component is present in 
those epochs close to the start of the outburst in the Balmer 
lines and a hint is present in the He llines. On 2007 April 
20 the dynamics of this component is similar to behaviour 
on April 15. By 2007 April 23 the emission component has 
faded into the absorption lines, which maintain their radial 
velocity semi-amplitude throughout. The systemic velocity 
of this narrow emission component is consistent with that of 
the He 11 emission, while the broad absorption components 
in the Balmer lines have a rather different 7 which is perhaps 
due to a slowly precessing disc. 
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Table 3. Velocity profile parameters on 2007 April 15 



Line 


7 


K 




Identification 




kms~^ it 3cr 


kms"^ ± 3(T 


±3a 




H/3 emission 


-16.6 ±0.8 


3.0 ± 1.6 


0.28 ±0.06 


SUngshot/WD? 


H/3 absorption 


22.5 ±4.4 


51.3 ± 6.5 


0.19 ±0.02 


Accretion disc 


H7 emission 










H7 absorption 


17.0 ±4.8 


43.6 ± 6.8 


0.18 ±0.02 


Accretion disc 


He II emission 


-19.4 ±4.7 


26.6 ± 6.8 


0.77 ±0.04 






0.0 0.5 1.0 1.5 2.0 

Phase 



Figure 4. Top panel shows the radial velocity curve of the He II 
emission. The radial velocity curves of the H/3 profile (bottom) are 
fitted with a low amplitude emission component and a high am- 
plitude absorption component. Both are from the phase-binned 
2007 April 15 data, showing two orbital cycles. 



GW Lib is one of the few CVs that have been followed 
through the outburst spectroscopically. For comparison, the 
spect ral evolution of SS Cyg (figure 5 in iHessman et al.l 
1 19841 ) showed broad, double-peaked emission lines that are 
gradually being overtaken by the rising continuum, followed 
by the formation of line wings in absorption. This quali- 
tative behaviour is expected as the accretion disc makes a 
rapid change from a low M optically thin configuration to 
a liigh A/ optically thick fiow at the onset of the outburst. 
Near maximum light, SS Cyg shows on top of the absorp- 
tion from the accretion disc, narrow emission cores which 
decay again after outburst. These are very similar to the 
spectral changes presented here for GW Lib, though no dy- 
namical properties are available to compare the components 
quantitatively. Unfortunately, the very low binary inclina- 
tion of GW Lib makes it difficult to study the accretion disc 
dynamics during the main outburst period. The absorption 
dominated disc lines are unsuitable for Doppler tomogra- 
phy and we therefore cannot search for disc asymmetries 
such as the tidal spirals seen in other DNe (e.g. IP Peg: 
ISteegs. Harlaftis fc Horn'31 19971 : U Gem: lGrootll200ll) . When 
disc emission returns, little structure can be detected in its 
marginally double-peaked lines. 

The stationary emission component in the spectra, visi- 
ble only during the first couple of days of the outburst, does 



not obviously fit in with the typical components expected 
to dominate the line emission in a mass transferring binary 
system. The low velocity suggests a location near the cen- 
tre of mass in the orbital plane or along the axis through 
the centre of mass, perpendicular to the orbital plane. GW 
Lib is not the first DN to show these features. Several DN 
systems have been reported to show low amplitude strong 
emission components in the Balmer lines. 

^tccghs ct al. (1996) report low velocity emission com- 
ponents in the Ha and He 11 4686A hues of the DNe IP Peg 
and SS Cyg and proposed slingshot prominences from the 
donor star as a possible origin for these features since their 
known system parameters rule out an origin near either the 
WD or the surface of the donor star. The rapidly co-rotating 
donor star may form prominences and such magnetic loops 
would be pulled towards the WD and could potentially find 
an equilibrium in between the LI point and the WD. As 
the prominence material is illuminated by the disc during 
outburst, it becomes visible as an emission source which is 
co-rotating with the binary orbit but located in the region 
near the centre of mass and thus displays very little radial 
velocity. If we follow their recipe, the observed H/3/He 11 
ratio is consistent with a gas temperature in the 10,000 - 
15,000 K range. 

In GW Lib, an alternative origin of the low- 
velocity emission may be near the surface of the WD, 
if its orbital velocity (Ki) is very small. Low state 
AM CVn systems often show such n arrow emission com - 
ponents (|Morales-Rueda et al.l |2003| . iRoelofs et all 120061 ). 
The low mass ratio deriv e d fro m the late superhump by 
iKato. Maehara fc MonardI (|2008l ') together with the low or- 
bital inclination would indeed imply a very low value for 
Ki in GW Lib. We will revisit these interpretations and the 
implications for the system parameters in section [5] 



4 GW LIB POST OUTBURST 

Binary evolution tells us that the majority of the current 
CV population should have evolved towards short orbital 
periods with the mass donor star depleted to a low mass 
degen e rate brown dwarf (e.g. iKn iggc 2006;:, iKolb fc Ritte3 
Il992l . iRappaport. Verbunt fc Josa .1983. ). Not only has 
this graveyard of CV s been elusive for many years (see 
iGansicke et"al] |2009| ). our sample of short period CVs 
with reliable binary parameters, other than the orbital 
period, is very small. This is mainly due to the fact 
that despite the low mass transfer rates, the donor stars 
are extremely faint and it is very difficult to perform 
traditional radial velocity studies using photospheric 
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absorption lines from the donor. Most studies are con- 



and constrain system parameters fe.e. table 2 in Marsh 


2001 


. ISouthworth et al. 


20071. Neilsen. Steeghs. & Vrtilek 


2008 


Mennickent. Unda-Sanzana & TaDDert 20061. 


Szkodv. Desai & Hoard! 


l2000h. The Ca ll triplet at 



more than make up for its relative weakness compared to 
the Balmer lines. It has a lower ionisation energy than 
either hydrogen and helium and is thus capable of being 
excited even by cool sources of radiation, and its thermal 
width and pressure broadening are much smaller than for 
hydrogen leading to sharper, mo re easily detected spikes of 
emission (|Marsh fc Dhillonlll99'it l. Studies have shown that 
Ca 11 is also more accessible than the Balmer lines for the 
study of velocity gradients and turbulence. This is because 
the thermal velocity broadening is smaller in general and 
the Keplerian shear broadening starts to domi nate at highe r 
inclinations compared to the hydrogen lines l|Hornelll995h . 
These studies show a promising, but so far neglected avenue 
for emission line studies of CVs. The average normalised 
spectrum of GW Lib in the I-band obtained on 2007 July 
24 is plotted in Figure [5] and shows a strong Ca ll triplet in 
emission. 



4.1 Emission from the secondary 

In Figure |S] (middle panels), we compare trailed spectro- 
grams of H/3 (left) and Ca ll 8662A (right) as obtained 3 
months after outburst. Whereas the Balmer line shows a 
shallow double-peaked profile and a rather blurry trail, the 
Ca II trail is much sharper even at a lower S/N. It shows a 
clear double-peaked profile from the accretion disc and an 
S-wave moving in between. No such S-wave appears to be 
present in the Balmer lines. 

To highlight this component, we averaged the profiles 
in a frame co-moving with the S-wave and show these in the 
top panels of Figure [S] The double-peaked profile from the 
accretion disc in this frame is averaged out, but the narrow 
S-wave is clearly visible on top on the broad profile in the 
Ca II case whereas its absence is noted in the Balmer profile. 

The S-wave is sharp and its amplitude is smaller than 
the velocity offset of the accretion disc peaks, suggesting 
it may originate from the secondary star rather than from 
the interaction point between the infalling stream and the 
accretion disc, the hotspot, since the latter would have a 
velocity equal to or larger than the outer disc edge. 

Assuming the emission does indeed arise from the sur- 
face of the donor star, the ephemeris can be determined as 
the crossing from blue-to-red of its velocity curve, which 
corresponds to inferior conjunction. This gives: 

HJD = 2454307.36867 + 0.05332S 

This is used throughout the paper to define the orbital phase 
although the accuracy of the orbital period is not enough 
to give the outburst data discussed in Section 3 a definite 
orbital phase since the accumulated uncertainty is too large 
(Figure [7]). 



^ From The Atomic Line List 
http : /www.pa.uky . edu/~peter/atoinic/ 
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The radial velocity semi-amplitude of the emission is 
determined using two different methods. The first method 
is the fitting of a Gaussian to the S-wave line component 
followed by a radial velocity curve fit. The second method is 
the localisation of the emission peak in the velocity-velocity 
plane using Doppler tomography. 

For the first method, the spectra are phase-binned and 
a triple Gaussian profile with as variables the offsets, the 
peak heights and the fwhms of the profiles is fitted where 
two Gaussians trace the disc features and the third is for the 
donor S-wave. Both the peak height and the fwhm show 
a chaotic variability with a maximum amplitude of 8 per 
cent from the average value and were therefore fixed to their 
mean values, leaving just the velocity offset as fit parameter. 
The offset velocity for the S-wave emission is plotted in Fig- 
ure [7] together with a weighted radial velocity curve fit. For 
improved accuracy in the Ca ll triplet, we performed a joint 
fit to the 3 triplet lines with a common velocity offset. We 
then find a semi-amplitude of Kq^ n = 82.2 ± 4.9 km s~^ 
for the S-wave in the Ca ll profiles, see also Table |4] 

To complement our radial velocity fits, Dopple r maps 
are c reated with a Maximum Entropy Method (MEM: lMarshl 
I2OOII ) for both the Balmer and the Ca 11 8662A lines, see 
bottom panels in Figure|S]for mem map of the time-resolved 
spectra on 2007 June 24. Given the fact that the lines from 
the Paschen series overlap with the Ca 11 triplet lines, we al- 
ways ensured that the Doppler reconstructions fitted to both 
lines simultaneously. The underlying Paschen line distribu- 
tion is rather featureless, and the sharp S-wave component 
is only present in the Ca 11 reconstructions. Even if such a 
Paschen contribution is not included, the Ca 11 tomograms 
show very similar features. 

The Balmer map (bottom left panel) shows a diffuse 
ring of emission from the accretion disc. Due to the sharper 
line profile, the Ca 11 map (bottom left panel) shows a much 
better defined accretion disc ring with a significant asymme- 
try at phases ~ 0.10, ~ 0.35 and ~ 0.70. Furthermore, the 
S-wave component manifests itself as a sharp emission spot 
peaking at a velocity of ~ 85 ±5 km s~^. Its location relative 
to the disc ring is as expected for donor emission and thus 
the Call S-wave is revealing the secondary in GW Lib for 
the first time. 



4.2 Emission from the accretion disc 

With the outburst activity decaying away, the disc contri- 
bution to the overall luminosity decreases and significant 
regions that are optically thin start to re-appear. These 
give rise to the characteristic double-peaked profile visi- 
ble in both the Balmer lines and the Ca 11 lines. A sin- 
gle broad Gaussian profile fit to this feature in the phase- 
binned spectra together with a weighted radial velocity fit 
for both the Ca 11 and Balmer profiles at multiple nights 
gives an average value for the amplitude of the disc of 
K^cc ~ 19.2 ± 5.3 km s~^, see Table |4l Similar results are 
obtained using multiple Gaussian line fits and diagnostic 
diagrams with double Gaussians. This amplitude is signifi- 
cantly different from the amplitude of disc related features 
in the April 2007 and July 2004 data. The phase shift of both 
Ca II and the Balmer disc coincide but are not in anti-phase 
with the donor star with respectively 4> = 0.68 ± 0.04 and 
(f) = 0.73±0.04. They cannot directly track the WD and thus 
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Table 4. Velocity profile parameters on 2007 June 25 



Line 


7 


K 




Identification 




kms~^ ± 3(7 


kms~^ ± 3(7 


±3<T 




Balmer emission 


-7.6 ±5.6 


19.2 ±5.3 


0.68 ± 0.04 


Accretion disc 


Ca II emission 


20.6 ± 3.3 


19.8 ±4.2 


0.73 ± 0.04 


Accretion disc 


Ca II emission 


-13.1 ±3.7 


82.2 ±4.9 


0.00 ±0.01 


Donor star 




200 400 - - 200 

km s~ km s~ 



Figure 6. Comparison of the dynamics of the Balmer lines (H/3, left column) versus Ca II (8662A, right column). 
The top row show average profiles in the reference frame of the donor star S-wave. The sharp donor component 
stands out in Ca ii but the HP profile is smooth. The middle row shows the phase-binned trailed spectrograms. 
The S-wave can be traced in between the disc peaks in Ca II only. The bottom row plots the Doppler tomograms 
showing the superior sensitivity and sharpness in the Ca II lines. The S-wave is mapped to a sharp spot consistent 
with donor star emission. Overlaying the Doppler tomograms are the positions of the WD, the Roche Lobe of the 
donor star and the streams for a model with q = 0.062. The inset in the Ca II Doppler map shows the Roche Lobe 
on the Donor star in more detail for both a. q = 0.062 and q = 0.23 solutions. 
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-6.0 ± 5.0 km s ^ for the Ca ll map. The fits for the 
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Figure 5. The average normaUsed I-band spectrum of GW Lib on 
July 24, 2007 showing the Paschen series and the double-peaked 
Ca II lines. 
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Figure 7. Top panel shows the radial velocity curve of the H/3 
emission. The bottom panel shows rotational velocity curve from 
the Ca II 8662A S-wave emission component. Both are from the 
phase-binned 2007 June 24 data, showing two orbital cycles. 



H/9 give identical values but with a larger error. The low Vx 
is encouraging from the point of view of linking the disc cen- 
tre with the WD at (0, -Kx). The implied A'disc = -Vy = 6 
km s^^ is smaller than the semi-amplitude derived using the 
Gaussian fits. On the other hand they are formally consis- 
tent when considering our error estimates, and both point 
to a small K\. 



4.3 Systemic velocity 

The binary systemic velocity (7) is a quantity that should 
obviously be the same in all spectral lines, refiecting the 
fixed radial velocity of the binary system with respect to us. 
For different lines at different epochs, the systemic velocity 
was determined using two methods. 

The first method determines the systemic velocity as 
derived from the radial velocity curve discussed previously, 
see Tables |3] and [l] For the Balmer lines associated with 
the accretion disc 7 — —26 ± 6.7 km s~^ before outburst, 
during outburst 7 — 22.5 ± 4.4 km s~^ and after outburst 
7 — —7.3 ± 3.6 km s^^. The narrow H/3 emission during 
outburst has 7 = — 16.6±0.8 km s~^. For our Ca 11 profiles, 
we found 7 — 20.6 ± 3.3 km s~^ for the disc emission and 
7 = —13.1 ± 3.7 km s~^ for the donor star emission line. 
The latter is the most reliable value for 7 as the donor star 
has to move with the systemic velocity. 

The alternative method uses the dependence of mem 
maps on the assumed systemic velocity. Reconstructed emis- 
sion spots in the map have optimal sharpness at the correct 
value for 7, while reconstructions at significantly different 
values will broaden the features and introduce possible arte- 
facts. Ranging 7 between -50 — 1-20 km s~^ in steps of 5 
km s~^ gives the best reconstruction for 7 = — 15±5 km s~^. 

Our various estimates therefore point towards a 7 of 
— 15 ± 5 km s~^, in ag r eemen t with the measurements by 
ISzkodv. Desai fc Hoard! (|2000l ). The only outliers are the 
Balmer absorption line during outburst with a systemic ve- 
locity fit of 22.5 ± 4.4 km s~^ and the Ca 11 emission asso- 
ciated with the accretion disc with 7 = 20.6 ± 3.3km s~^. 
However, those fits are distorted as the resolution is not suf- 
ficient to properly resolve the donor star S-wave and the 
disc emission peaks near phases 0.25 and 0.75. Subtraction 
of the S-wave changes and fitting the disc emission with a 
two-Gaussian profile with an fixed offset between the Gaus- 
sians gives a lower value of 7 = 7 ± 2 km s^^. 



the radial velocity amplitude is disturbed by some residual 
disc asymmetry and cannot be straightforwardly connected 
to Kx. 

The accretion disc in the Doppler maps of H/3 appears 
rather featureless with little azimuthal structure, while the 
Ca II trail and mem map show lower intensity in the disc 
around phases ~ 0.10, ~ 0.35 and ~ 0.70. The centre of 
emission of the discs for both maps was determined by look- 
ing at their centre of symmetry. To determine the optimal 
centre of symmetry of the disc, we subtracted a symmet- 
rical component from both the Ca 11 and the H/3 Doppler 
map centred at T4, Vy in the ranges of -100 — 1-50 km s^^ 
making grid steps of 2 km s~^ and insp ected the residuals 
in th e lower half of the map (see also ISteeghs fc Casarei 
|2002| '1. The best fits are found for K = 0.0 ±5.0 km s"^ and 



5 SYSTEM PARAMETERS 

In the previous section we have seen the superior sensitiv- 
ity of the Ca 11 emission lines over the commonly observed 
Balmer lines. Not only do they provide a much sharper view 
on the accretion disc emission, the key result was the pres- 
ence of an emission spot from the donor star. This has pro- 
vided the first ever proxy for its orbital velocity. 

Post-outburst radial velocity curves together with 
Doppler maps have given us estimates of the semi- 
amplitudes for the different components in the binary sys- 
tem. From both the Balmer and Ca 11 double-peaks, the 
radial velocity of the accretion disc was found to be K^^c = 
19.2 ± 5.3 km s~^. However, this value may be biased by 
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disc asymmetries and the donor star S-wave. The observed 
phase-shift between the disc peaks and the donor does in- 
deed suggest such a bias to be present, while the centre of 
symmetry searches returned lower Ki estimates. We there- 
fore consider -ftTacc from the disc peaks as an upper limit on 
the radial velocity of the WD: Ki < A'acc- 

The donor star is detected in the Ca ll emission as a 
third peak on top of the accretion disc emission and gives 
Kam ~ 82.2 ± 4.9 km s~^. Considering the faintness of the 
low mass donor star, this is unlikely to be powered by chro- 
mospheric emission. In any case, the Ca ll triplet in chromo- 
spherically active stars is generally seen in absorption (e.g. 
iKafka fc Honevcuttl200^ ) and thus cannot explain the emis- 
sion seen in GW Lib. A more likely origin is photo-ionisation 
in the irradiated hemisphere of the donor star facing the 
WD. The measured radial velocity semi-amplitude is then 
an underestimate of the true orbital velocity of the donor 
star {K2 > Kem) since not all of the Roche lobe contributes. 
We calculated the magnitude of this bias, referred to as the 
K-correction, by simulating the expected emission profiles 
from an irradiated Roche lobe filling star. Synthetic pro- 
files were calculated for relevant binary parameters includ- 
ing possible shielding of the equatorial regions of the donor 
by a vertically extended accretion disc. 

The peak intensity of the donor star S-wave does not 
vary significantly as a function of the orbital phase, and in- 
stead scatters around a mean value with an RMS of 8 per 
cent, which is consistent with photon noise. The model pro- 
files discussed above were also used to generate expected 
orbital light curves of the donor star emission as a function 
of binary parameters. Model light curves, assuming that the 
emission is optically thick and thus should be weighted with 
the projected area, produce orbital modulations above what 
is observed. However, if we consider more isotropic emis- 
sion, the modulation disappears at low orbital inclinations. 
This optically thin model is more appropriate for emission 
stimulated by photo-ionisation and thus the observed lack 
of variability is consistent with an irradiated secondary ob- 
served at low inclinations. 

We can now derive conservative limits on the radial ve- 
locities of the stellar components in the system using vari- 
ous methods and spectral lines. They can be summarised as 
K2 > 82.2 ± 4.9 km s"^ and Ki < 19.2 ± 5.3 km s'^. These 
limits translate into a hard upper limit for the mass ratio of 
q < 0.23 regardless of the magnitude of the K-correction or 
the Ki overestimate. 

After the April 2007 ou tburst, two detections of su per- 
humps have been reported. ICopperwheat et al.1 (|2009l ) de- 
tect a periodicity with a period excess of 4.12 min which im- 
plies a mass ratio of g ~ 0.211 if interpreted as a superhump 
and the empirical supe rhump-excess mass ratio relation from 
[Patterson et al.l ()2005l ) is applied. This is close to our hard 
upper limit for the mass ratio. However, the authors them- 
selves question the character of the humps and the implied 
donor stars mass in combinatio n with the WD mass from 
seismo logy. On the other hand, iKato. Maehara fc MonardI 
(|2008l l report a superhump period of 0.053925(4) days and 
a late superhump period of 0.054156(1) days which imply 
an extreme mass-ratio of g = 0.062. 

We looked at a simple model for the relation between 
M2 , Ki and K2 as function of q for a range of Mi . Combin- 
ing this with the above upper and lower limits for the radial 



Table 5. Deriv e d sys tem parameters for GW Lib using Mi from 
iTownsIev et al.l |2004i). The left hand column lists the binary pa- 
rameters derived from the formal solution based on the measured 
limits for Ki and K2- The right hand column lists the param- 
eters when combining the mass ratio inferred from superhumps 
with the calculated K2 from the Ca II Kcm- 





q < 0.23 


q = 0.062 


Ml (Mq) 


1.02 


1.02 


K2 (kms-i) 


> 82.2 


100.8 ±7.1 


Ki (kms-i) 


< 19.2 


6.2 ±0.4 


M2(Mq) 


< 0.23 


0.063 


vi (kms"-*-) 


< 122 


34 


V2 (kms~^) 


> 488 


548 


in 


> 9.4 


10.6 ± 0.8 


SpT 


M3.9-M5 


M8-T 



velocities of respectively the WD and donor star rules out 
systems with Mi < 0.75A'Iq. As the system is tidally locked, 
the rotation period of the secondary is equal to the rotation 
period of the system. Emission lines are broadened by the 
rotation of the secondary star where the surface velocity can 
be approximated by looking at the co-rotating Roche Lobe 
surface. Comparing the measured fwhm of the emission line 
to the theoretical value as a function of the mass ratio can 
provide another limit on q and Mi. The observed fwhm of 
Ca II line emission peak in each bin ranges between 1.0 and 
2.2A due to contamination with the double-peaked profile. 
The FWHM is very close to the resolution of the data, 0.617A, 
which equates to a n sin i of 21 km s~^. The profile is thus at 
best marginally resolved. From the upper limits for the mass 
ratio and the constraints given by the measured emission 
line broadening we get a rough estimate of the WD mass at 
l.OMp) ±0.25, in agree ment with the mass determination by 
iTownslev et al.l l|2004h . Note that their window of solutions 
is constrained by the UV-fiux of the WD in combination 
with the measured dista nce to GW Lib. For a parallax dis- 
tance of IO4I20PC (from Frhorstensenl |2003| ) they calculate 
M l = 1.03 — 1.36M(;), supe rseding the previous esti mates 
by 'Szkod v et all (120021 ') and iThorstensen et al.l l|2002l ') that 
were based on a larger distance. For further analysis we will 
use the asteroseismologicaly suggested mass for Mi. 

When con sidering the empirical CV donor sequence of 
iKniggel l|2006l ) at the orbital period of GW Lib, we find a 
typical pre/post period bounce secondary star mass of AI2 — 
O.O64±O.OO1M0 or M2 = 0.060±0.001Mq , respectively. The 
small difference is due to the orbital period being very close 
to the bounce period. Knigge also calculates the expected 
absolute magnitudes for GW Lib in the infrared bands using 
its parallax distance, which implies a low mass secondary 
star of O.O8O±O.OO5M0. 

5.1 Formal solution: q < 0.23 

The formal scenario assumes that corrections that need to 
be applied to our Ki and K2 constraints are modest. Thus 
implying that the mass ratio is near its limit: q = K1/K2 < 
0.23. The right hand column in Table [5] list the derived sys- 
tem parameters for this ratio. IThorstensen et al] (|2002l ) cal- 
culated the inclination based on the fwhm of the shallow 
peaked Balmer emission lines in quiescence compared to the 
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FWHM of WZ Sge and found i ~ 11°. This is in agreement 
with our calculations in combination with the nominal solu- 
tion for Ml = 1.02 M© and confirms that GW Lib is indeed 
observed at very low inclination. 

For q < 0.23, the correction for K2, as derived from our 
models, is at most AisTmax ~ 35 km s^^ giving a maximum 
value for K2 of 118 ± 5 km s~^. 

The main problem with a high value for q near our limit 
is that the secondary would then need to be a M2 ^ 0.23 Mp -, 
mass star with a spectral type around M4-M5 (|Coxll2000l ). 
In that case, spectral features should be visible in a high 
resolution, high signal to noise red spectrum. We cannot 
identify any of such in our post-outburst Lband spectra nor 
in the pre-outburst Magellan data. Additional issues are the 
observed IR magnitudes of GW Lib which imply a faint 
donor and a main-sequence object near 0.23 Mq would not 
fit within the Roche-lobe. All adding to the conclusion that 
the mass-ratio in GW Lib is likely not near this maximum 
value. 



5.2 Superhump solution: q = 0.062 

Measuring the late su perhu mps in GW Lib, 
iKato. Maehara fc MonardI (|2008h found an extreme 
mass ratio for the system of g = 0.062. Which, when com- 
bined with the WD mass fr om seismology and the empirical 
donor star mass from the iKniggd (|2006l ) sequence, proves 
an attractive solution within the boundaries of the formal 
scenario. We calculated the magnitude of the K-correction 
at this mass ratio using our irradiated secondary model. In 
order to produce the observed Kcm ~ 82.2 ± 4.9 km s~^, 
we find a true K2 — 100.8 ± 7.1 km s"'^ where the error 
is calculated by propagating the error on Kem- This 
error dominates over effects such as the assumed level of 
equatorial shielding or the optical thickness of the emission. 
The implied system parameters based on this assumed q 
and calculated K2 can be found in the left hand column of 
Table E] 

The small mass ra tio scenari o would correspond to a low 
meiss donor with a M8 (|Coxll2000l ) or T spectral type l|Kniggel 
I2OO6I ) if the donor is close to its main-sequence configuration. 
This would agree with the infrared magnitude constraints 
mentioned previously and fits perfectly with GW Lib being 
close to the period minimum. 

Interestingly, the derived Ki in this scenario is close to 
the measured radial velocity of the narrow Balmer emission 
component in outburst {K = 3.0 ± 1.6 km s~^) suggesting 
that this component may trace the WD when irradiated by 
the accretion disc. Similar sharp emis sion components from 
the WD arc seen in AM CVn systems (jMorales-Rueda et al.l 
[2003, Roclofs et al..,2006i ). However, a component arising at 
the surface of a WD is expected to be gravitationally red- 
shifted but the systemic velocity for this lines component 
was —16.6 ± 0.8 km s~^ which is similar to the systemic 
velocity of the system and cannot be biased by the Stark 
effect as is the case in AM CVn systems. Unfortunately the 
emission is only a transient feature and the constraint on 
its K will be difficult to improve upon. The origin of the 
component and its possible connection with the WD thus 
remains unclear. 



6 DISCUSSION 

We presented time-resolved optical spectroscopy of GW lib 
during a large number of epochs spanning before, through- 
out and after the 2007 outburst. We studied the long-term 
evolution of the spectral features tracking large changes in 
the accretion geometry and intensity. 

Pre-outburst spectroscopy obtained in 2004 data show 
clear features of the accreting WD including a detection of 
Mg 11 absorption. Here, the low inclination of GW Lib has 
the advantage over high inclination system as it gives narrow 
disc emission lines instead of wide profiles, potentially re- 
solving the line from the nearby He I line. Higher resolution, 
high signal to noise spectra around this line in quiescence 
could provide a dynamical trace of the accreting WD and 
give a directly measured mass constraint to complement the 
constraints indirectly obtained from seismology. 

During the outburst, we initially see the optically thick 
accretion disc dominate through broad absorption features. 
Their radial velocity does not trace the WD and show semi- 
amplitudes of order 50 km/s. A peculiar sharp emission line 
component was found in the Balmer lines that is effectively 
stationary. This could be some low velocity outflow in the 
z-direction or a component near the WD, but its mean veloc- 
ity is similar to 7 whereas a layer near the high gravity WD 
should be gravitationally red-shifted. As the system fades 
back towards quiescence, accretion powered double-peaked 
emission profiles appear. The absorption associated with the 
accretion disc continues to weaken until the system returns 
to a semi-quiescent state with strong double-peaked emis- 
sion lines. 

The Ca 11 triplet in GW Lib shows that these much 
neglected lines of CV spectra could be an interesting win- 
dow to search for signs of both donor stars and accretion 
disc structures even in cases where the Balmer lines show 
no signs of the donor star whatsoever and the disc itself is 
barely resolved. GW Lib is a particularly challenging object 
in that sense due its very low inclination and thus small 
projected velocities. The Ca 11 triplet in emission could be 
resolved into several components and the secondary star was 
discovered in emission moving in between the sharp double- 
peaked emission from the accretion disc. Doppler tomogra- 
phy and radial velocity profiles fits of the Ca 11 lines pro- 
vide a semi- amplitude of Kcm = 82.2 ± 4.9 km s~^ for this 
donor star feature and indicates a disc centre of symmetry 
at iodise = 6 ± 5 km s^^. The disc is also visible in the H/3 
maps but with less detail and sharpness and no donor star 
contribution is seen. Contrasting these lines highlights the 
diagnostic advantages provided by the Ca 11 triplet. 

Based on previous studies together with the limits on 
Ki and K2 provided in this work, the allowed range of binary 
parameters were explored. While our dynamical limits place 
a hard upper limit on the binary mass ratio of g < 0.23, 
we favour a significantly lower value. A mass ratio near q ~ 
0.06 is in accordance with estimates based on the detected 
super-hump modulations, the constraints on the faintness, 
and thus mass, of the donor star and the indications that 
Ki is very small. Given such a mass ratio, our measured 
Kcm implies K2 = 100.8 ± 7.1 km s~^ for the donor star 
component when applying the relevant K-correction. The 
implied Ki — 6.2 ± 0.4 km s~^ is then also close to our 
measured disc centre of symmetry. The combination of a 
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WD mass near the value suggested by the pulsations and a 
low mass donor near the empirical sequence of an evolved 
CV near the period bounce appears to be consistent with 
the observational constraints to date. 

Whether the spectral features observed several months 
after the tail end of the 2007 outburst are persistent remains 
to be established. As the WD continues to cool, donor star 
irradiation may be less effective in exciting the strong Ca ii 
lines we observed. Nonetheless, high resolution spectroscopy 
resolving these lines at good S/N levels appears to be a vi- 
able tool to expand our knowledge of the binary parameters 
of short period cataclysmic variables and their faint donor 
stars. 
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